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Anisotropic neutrino emission during the neutron star formation can be the origin of the observed 
proper motions of pulsars. We derive a general expression for the momentum asymmetry in terms 
of the neutrino energy flux gradient, and show that a nonvanishing effect is induced at the lowest 
order by a deformed neutrinosphere. In particular, this result is valid for a neutrino flux transported 
through a spherical atmosphere with constant luminosity. 
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^ • I. INTRODUCTION 

O 

10 ' Observations show that pulsars have peculiar proper motions. They have very high translational velocities with 
OA , respect to the surrounding stars, with a mean value of 450 km/ s and up to a maximum of about 1000 km/ s [|l|. This 
suggests that some kind of impulse (kick) happens during the birth of the neutron star. Different mechanisms have 
been proposed to explain the kick, but most of them have difficulties to produce the large observed velocities. 

Neutrinos carry away almost all the energy released in the gravitational collapse (« 3 x 10^'^ erg), taking with them 
' a momentum ~ 100 times the momentum associated with the spatial motion of pulsars. Therefore, a 1% anisotropy 
in the momentum distribution of the outgoing neutrinos would suffice to account for the translational kick. 
\^ An interesting mechanism to explain the asymmetric neutrino emission from a cooling protoneutron star has been 

proposed by Kusenko and Segre ||]. It is based on the matter neutrino oscillations in the presence of an intense 
magnetic field. The emission surface of the electron neutrino is located at a radius larger than the one corresponding 
to the muon or tau neutrino. Under suitable conditions, a resonant transformation v,, Vt can take place in the 
region between the boundaries of the electron and the tau neutrinospheres. The v^, are trapped by the medium, but 
C^' the Vt produced in this way are outside their neutrinosphere and free to escape from the protostar. Consequently, 
Q the surface of resonance acts as an effective tau neutrinosphere. If there is a magnetic field, or another non-isotropic 
?H ' effect, this surface of resonance becomes distorted and an anisotropy in the energy fiux is generated, causing a kick 
c/3 , to the protostar. 

' Doubts about the effectiveness of the above mechanism have been raised by Janka and Raffclt ||]. According 
^ \ to them, no effect is generated at lowest order because it is not justified to calculate the flux asymmetry from the 
temperature variation around the surface of resonance. The neutrino luminosity in the protoneutron star is controlled 
by the core emission and is not affected by local processes in the atmosphere, where the flavor transformation occurs. 
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^ In support of their argument, Janka and Raffelt use the Eddington model for a plane-parallel stellar atmosphere jj] 
■ ■ ■ ' to estimate a residual asymmetry that is induced by higher-order corrections. 

In this work we reconsider the problem of the neutrino oscillation mechanism for pulsar kicks. In fact, a neutron 
protostar is a very complex system, mainly formed by interacting nucleons, electrons, photons and neutrinos, rotating 
quickly and with a strong magnetic fleld. However, for our analysis a simplified and perturbative description is 
suitable. The star can be considered as constituted by two gases, the nucleon gas and the neutrino gas. From the 
hydrodynamic point of view electrons and photons can be ignored. The dynamics of these gases are controlled by the 
gravitational field, the magnetic field, the neutrino-nucleon interactions and the hydrostatic equilibrium equations. 
At the lowest order in this perturbative approach, and neglecting the star rotation, the dynamics of the nucleon gas 
is dominated by the gravitational field and the equilibrium hydrodynamic equations, and it is thus described by a 
spherical distribution. Besides this we have the neutrino gas, interacting with the nucleon gas and the magnetic field. 
This last interaction acts in particular on the neutrino mass columns, breaking their isotropy. This breaking alters 
the hydrostatic equilibrium equations, producing a non-isotropic distortion in the nuclear matter. In the following we 
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neglect these contributions, but even so there is a non vanishing kick effect induced by the geometrical deformation 
of the resonance surface, provided that it acts as an effective emission surface. 

We parametrize the resonance surface as R ~ Rr + 6 cos 6, and develop the expressions up to first order terms in 
S/Rr- From this approach we derive an expression for the fractional momentum asymmetry in terms of the spatial 
derivative of the energy flux. Using this result, we show that the distortion of the resonance surface by the presence of 
a magnetic field generates a geometrical asymmetry in the neutrino emission, even in the case of a constant luminosity. 
To illustrate this effect, we consider two simple self-consistent models for a spherical protostar atmosphere, which 
satisfy the energy flux conservation. In particular, in one of the models we use the Eddington approximation adapted 
to the spherical geometry. The required magnetic field results to be one order of magnitude higher than the value 
derived in the original articles . 

To make the comparison with the existing literature easier, we restrict the discussion to the standard mechanism 
of neutrino oscillations between massive active neutrinos, and a magnetic-field induced deformation of the resonance 
surface. Nevertheless, our approach can be straightforwardly extended to other situations [^, and in particular to 
oscillations produced by a violation of the equivalence principle (VEP) In the case of VEP no magnetic field 

is needed to deform the resonance surface 

In the following section we examine the effect of the deformation of the resonance surface on the neutrino energy 
flux, emphasizing the relevance of the geometrical variation of flux to produce an asymmetric momentum emission. 
In Sections III and IV, we apply the results of Section II to the Eddington and the polytrope neutrinosphere models 
and estimate the magnitude of the magnetic field required to explain the observations. The last section presents some 
general conclusions. 



II. SURFACE OF RESONANCE AND NEUTRINO ENERGY FLUX 



We consider oscillations between two neutrino flavors, say Ve and i^r, in the interior of a protostar. Neutrinos have 
an average energy E = k = |k|, which depends on the radial coordinate r. In the absence of a magnetic field, or any 
isotropy- breaking interaction, the resonant transformation takes place on the surface of a sphere of radius R^, given 
by the condition 

^^COS20 ^ V2GFNer , (1) 

where — k[Rr) and N^.^. — Ne{Rr)- Here, 9 is the vacuum mixing angle, Gp is the Fermi constant, and Am^ = 
m| — rn^ is the difference between the square mass of the neutrinos. We assume that the number density of electrons 
iVe(r) is proportional to the baryon density p{r), N^ir) ~ ^pp(r). 

In the absence of a magnetic field the emission surface is spherical, the outgoing energy flux is radial and there is 
no kick, but the presence of a magnetic field, described here by a uniform field B, distorts the surface of resonance. 
In the simplest case this is now defined by a function R{'d) ~ Rr + (5 cos??, where 5 < Rr and d is the angle between 
the vector position of a point at the surface and the direction of the magnetic field of the protostar. The distortion of 
the resonance surface leads to a modification in the outgoing energy flux. The kick is characterized by the fractional 
momentum asymmetry factor 
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where Fg is the outgoing electron neutrino energy flux at the element of area da of the emission surface. The integrals 
in the denominator and numerator give the total momentum lost by the protostar per unit of time and its component 
along the direction of the magnetic field, respectively. The factor 1/6 comes from the fact that only the electron 
neutrino contributes to the energy flux asymmetry. 

To compute the asymmetry factor it is necessary to take into account the structure of the flux at the resonance 
surface, considered as an effective emission surface. The details of this analysis are presented in the appendix. The 
basic information is given by the distribution function for neutrinos in thermal equilibrium with the medium and 
satisfying the diffusion approximation 

^ - 1 - df^"^ 
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where k = kCl is the momentum of the neutrinos, f^'' = (l + e^'^"^"^/-^) is the neutrino distribution function 
at equiHbrium, is the chemical potential, and the factor A comes from the differential cross section for neutrino 
reactions. The energy flux is defined by this distribution function according to 



F.(r) 



(2^ 



k/.(r,k) 



(4) 



At the interior of the resonance surface the energy flux is radial (F = Ff) and results from the diffusive part of 
the neutrino distribution. As discussed in Ref . Q , in a neutron protostar the neutrino luminosity Lc is governed by 
the energy loss from the core. Throughout the neutrinosphere, within the resonance surface, the luminosity can be 
assumed as not being dependent on the radial coordinate, and thus as satisfying F(r) — Thus, below the limit 

surface, the flux is purely diffusive, and is given by the constant luminosity condition. Once the energy flux above the 
limit surface is computed using the expression (0), we find that it has a radial component, associated to the diffusive 
part of the neutrino distribution, and a normal component, due to the isotropic part of the neutrino distribution 



Ff, n + Ff f , 



(5) 



where fi and f are unit vectors along the normal and radial directions respectively, at the considered point on the 
resonance surface. The local flux conservation (V • F = 0) implies that, at the resonance surface, the normal outgoing 
flux is equal to the diffusive outgoing flux, and that both are equal to one half of the diffusive flux calculated just 
below the resonance surface 



Ff{'9) ^ F,,(d) = \F,(Rr 



Scosi}) ~ ^F^{Rr) (1 
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where h. 



1 dF 



In the integrals of Eq. (^), da is the element of area on the distorted surface of resonance and does not coincide 
with the element of area dur on the sphere of radius Rr. These areas are related as follows: 



da — 



1 dR 
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In addition, we have f • B = cosz9, and 
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From Eqs. (|^)-(|^), to first order, we observe that only Fn, the component of the energy flux normal to the 
resonance surface, gives a non null contribution to the integral in the numerator of the asymmetry factor, and it is 
thus responsible for the kick. It results 



k 36 ^ 



(9) 



Taking into account the energy emitted by the protostar in form of neutrinos, this ratio must have a value of the 
order of 10~^ to produce the observed kicks. This expression clearly shows that the existence of a kick requires a 
nonvanishing gradient of the flux, i.e. hp^ ^ 0. Using hp^ = — 2_R^^, we obtain 
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(10) 



This means that S must be of the order of Rr/6 to produce the required kick. 

To calculate S, let us remember that the index of refraction of the neutrinos is modified by the presence of an 
external magnetic field |lO[ . This fact affects the flavor transformations of mixed neutrinos producing an anisotropic 
contribution to the resonance condition |pl]] . For neutrinos propagating through a degenerate electron gas, the 
resonance condition becomes 
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where e is the electron charge, B = |B|, and k and A^e are evaluated at R{'d). In Eq. (11) it has been assumed 
that the weak-field limit is satisfied, i.e. 2eB <C (STr^A^e)^^'^- More general features of the neutrino propagation in 
magnetized media, incorporating the effect of strong magnetic fields, have been considered by several authors ||l^. In 
our case, the weak-field condition means B ^ 5 x lO^^G and, in fact, is marginally satisfied in the models discussed 
below. However, taking into account the ambiguities in the values considered for the parameters of a protoneuton 
star, a more accurate computation is not necessary. This could be meaningful in the context of a more precise and 
detailed model, which is beyond the scope of this article. 



By writing k = kr + Sk and 



SNe, from Eqs. (0) and ^ ) we obtain the relation 



Sk cos 261 ~ -V2Gf SN^ 



71" 



3jVe 



1/3 



B cos-d 



(12) 



There are two contributions to SN^. One is due to the geometrical distortion of the surface of resonance, assuming 
that the different profiles of the neutrinosphere remain unchanged. The other comes from the distortion of the profiles 
of temperature, pressure, density, etc., induced by the geometrical distortion. The equations that define the model 
give the relation between the last higher-order contribution and the deformation of the surface of resonance. In this 
work we will only consider the first order contribution. Thus, we have 



dr 



5R = hj^^NerS cos l} 



Analogously, we also have 



6k = 



dk 
dr 



6R = hf^ ^krS cos d 



Inserting Eq. (|l|) and (|T|) into Eq. (|l|), we get 



e / 3 
2 W 



1/3 



N-^/^B- 



(13) 



(14) 



(15) 



If we assume that the electron neutrinos are in thermal equilibrium with the stellar medium, the average energy of 
the emitted neutrinos is proportional to the temperature at the emission point, k = j^^j^T ~ 3.15T. In such a case, 
from Eq. (P we get 

Am^cos2e~^^^ . (16) 



To have the resonance within the electron neutrinosphere, PrT^ must be larger than the corresponding value at the 
surface of the neutrinosphere, p^^T^^. For an ideal gas this simply means that the pressure at the resonance must be 
larger than the pressure at the surface of the neutrinosphere. The z/e trapping density is p ^ lO^'^ g cm~^, and hence 



Am^ cos 29 > 



4.3 X 10"^ X 



(17) 



where the temperature Ty^ is given in MeV. For T^^ ~ (3 — 5) MeV, it requires m^^ > 100 eV. Now, hk = hr and, 
together with h^^ = hp, we obtain 



3eB 
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hj, + hp 



(18) 



In general, to compute 6 and hp (or Rr in the case of Eq. (pO|)) a model for the neutrino atmosphere of a neutron 
protostar must be specified. This will be done in the next sections, where we examine two analytical models that 
are meaningful up to the neutrinosphere, where the neutrino transport equation holds. In both models the neutrino 
transport is in the diffusion regime and the luminosity is independent of r, but they differ in the assumed properties 
for the medium. In the first one, the Eddington model, the medium is an ideal gas of nucleons, while in the second 
one it is a polytrope gas. 
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III. THE SPHERICAL EDDINGTON MODEL 



The Eddington model gives a simple and physically reasonable description of a neutrino atmosphere, locally ho- 
mogeneous and isotropic. For a plane geometry the model was developed by Schinder and Shapiro and here we 
extend it to the spherical geometry. 

For neutrinos and antineutrinos in thermal equilibrium with the medium, and satisfying a transport regime con- 
sistent with the diffusion approximation, the energy density, the energy flux, and the stress tensor of neutrinos with 
momentum k are [Q: 



27r2 



1 fc3 



eg 



1 



(19) 
(20) 
(21) 



We assume that A = npk , with k = 5.6 x 10^ erg~ cm . In the interior of the neutrinosphere we can consider 
that there are two perfect fluids. One of them is constituted by nonrelativistic nucleons of mass m„, with a density 
p, and the other by ultrarelativistic neutrinos and antineutrinos, with a vanishing chemical potential, = = 0. 
Photons and electrons are of course present and, in fact, electrons make the relevant contribution to the effective 
potential in the case of matter oscillations between active neutrinos. However, we can ignore both of them for the 
hydrodynamic description of the system. Therefore, the relationships between pressure P , density of energy e, and 
temperature T are 



P = 



m 



24 



P 



7a 



rp4 



(22) 
(23) 



where a ~ 2.09 x 10'*^ erg ^cm ^ is the Stcfan-Boltzmann constant. Taking into account the gravitational field (j) of 
the star, we have for P 

VP = -{P + p)\7(l3 . (24) 
Putting all these relations together, we obtain the set of equations which describes an isotropic neutrinosphere 

(25) 



U{r) = \aT\T) , 

36 Kp(r) dr 



T\r) 



P{r) 



^T(r) + |^T4(r) , 
m„ 24 



P{r) + p{r) +^-^T\r) 



GM{r) 



(26) 
(27) 
(28) 



with M(r) = 47r df f"^ p^(f) , where p^ is the total density of mass. 

In the region where a resonant transformation could happen the baryon density is p ~ (10^^ — 10^^) g cm^^ and 
T ~ (3 - 10) MeV. Thus, we have p^^ = ~ (lO^^ - lO^^) p and P ~ (lO^^ - lO^^) p. Therefore, Eqs. ^ 

and (Ba) reduce to 



P = 

dP 

dr 



pT 



-P- 



GM{r) 



(29) 
(30) 



These last equations together with Eq. ( pq ) describe an ideal gas of nucleons at hydrostatic equilibrium, with an 
energy flux given by the transport of neutrinos. 
Prom Eqs. (p|) and (|o|) we have 



5 



where is the protostar mass enclosed by the resonance sphere. Therefore, from Eqs. (|l|), ®, and (|l|), the 
necessary magnetic field results 



B ~ ^ " . (32) 

At this point, assuming that the resonance occurs near the surface of the neutrinosphere, i.e. Rr — Ru^, we can make 
an estimation of the order of magnitude of B for a typical protoneutron star. Then, taking Rr ~ 30 fcm, Mr — Mq, 
Tr ^'a MeV, and pr — 10^^ g cm~^, to have Ak/k ~ 0.01 the magnetic field must be i? ~ 10^® G . Of course, these 
parameters are not independent and a more careful discussion of the model is necessary. 

The neutrinosphere is defined by four functions: pressure P{r), temperature T{r), baryonic density p{r), and 
energy flux F(r). Up to this point we have only three independent equations relating these functions. A complete 
specification of the system requires a fourth relation. A simple analytical model that satisfies the requirement of a 
constant luminosity is the Eddington atmosphere 0|. This model is defined by Eqs. (^), (p9|), and (|30|), plus the 
hypothesis of the energy flux conservation \/.F = 0. For an isotropic fiux, the additional assumption leads to: 



^(^^^)-0, (33) 



dr 



which simply means 



475-7-2 ' 



(34) 



where is the luminosity of the protostar. 

From Eqs. ( |30| ) and (^4|) we have for the mass density 



2tt ^dT ^^^^ 
9 kLc dr ' 

and replacing this expression in Eq. ( |30| ) we arrive at 

dP 2^ G dT 

— = —M{r)T—. (36) 

dr 9 kLc dr 

To find the solution of the structure equations for this model we use the following procedure. First we define a 
reduced effective mass m(r), given by 

dT r dT 

dr M{r)T — = MMr) / dr T — , (37) 
dr Jr^ dr 

where Mc is the mass of the core. In terms of m(r) the solution of Eq. ( ^ ) is immediate. Defining etc = ^P^-^""^"'^'' 
we have 



9 kLcPc 



where 



a{r) = 1 5— . (39) 



ac'm{r) 

From Eqs. (|29|) and (^8|) we can express the density in terms of the temperature as follows 



pM = Ar^f^-a) , (40) 



1 - a r V^c 

and using this result in Eq. (p3) we obtain a first order differential equation for T 
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^■^4^- According to this equation the slope of the temperature at Rc is independent of the function 



where Ac = 

a(r), = — . If we refer to an ideahzed neutrinosphere, where this model would apply in the whole space, the 
interesting solutions correspond to infinite protostars where the temperature has an asymptotic behavior for r ^ Rc^ 
such that the temperature tends to Ts ~ \faTc. Thus, the function a(r) varies in the range 1 — a^^ < a < {Tg/Tc) 
for Rc < r < oo. 

This system of equations has no analytical solution when a is a function of r, and in general there is no perturbative 
expansion that gives a good approximate solution at every point within the neutrinosphere. To find an approximate 
solution let us consider the differential equation ( pj ) with a constant. In this case, for Tc > T > an analytical 

(implicit) solution is given by 



ln(^-l^')-ln 



1 - 



1 



1 



^-1 
a \ r 



(42) 



If we replace the constant a for a (well behaved) function of r, then the above expression still satisfies Eq. ( [4l| ) at 
r — Rc- For an infinite atmosphere, a good approximation to the exact solution is given by Eq. (^), with a now a 
function of T(r): 



1 T 
a(T) = l ^ l-TT 



(43) 



where A ■ 



and Ts denotes the temperature at the surface of the electron neutrinosphere, which 



is assumed to lie in the asymptotic region. From Eq. (p3|), we see that a{Ts) = (Tg/Tc) and a{Tc) = 1 — aj^, and 
thus it fits the extreme values of a{r). 

The surface of the electron neutrinosphere corresponds to a density ps ~ p^^ . Assuming that the temperature at 
this surface is close to Ts, from Eq. (M) we obtain 



T. 1 



Tr 



2Ts - ATc Pc 



(44) 



with p^^ <C Pc- The radius of the neutrinosphere is obtained by evaluating Eq. ( p2| ) in T^^. 

To illustrate the predictions of the model, we adopt a neutron protostar with reasonable values for its parameters. 
The core is assumed to have Mc = Mq = 1.13 x 10^" MeV, Rc = 10 km, Lc = 9.5 x lO^^ erg s-\ pc = lO^"* 



g cm 



and Tc = 40 MeV. We take the surface of resonance as defined by a density Pr — W^^g cm~^, while a numerical 
estimation for the asymptotic value of the solution of Eq. (^ij) gives Tg = 4.8 MeV. Inserting these values into Eqs. 
( ^ ) and (^) we get R„^ ~ 2.7Rc, and a numerical computation of the total mass of the star from Eq. (|4^) gives 
Ms — 1.4 M(7) . Then, for a resonance region lying near the surface of the electron neutrinosphere, from Eqs. (|4T 
(EO), and dl) we have 



-Ac 



T^RcPr 



.04- 



=Rr- 



TrPc 
TcPr 



A- 



Tr 



Rr 



19- 



Rr 



(45) 
(46) 



According to Eqs. ( p^ ) and (|To|), we see that a Ak/k of order 0.01 can be obtained with i? ~ 3 x lO'^^ G. This value 
is in agreement with the estimation done by means of Eq. (|3^) . In general, for more extended and hotter Eddington 
protostars smaller magnetic fields are needed, as can also be seen from Eq. (^. 

The strength of the magnetic field we have obtained is somewhat higher than those estimated in previous works on 
the subject |^JHl. However, it is important to note that for us B is at least an order of magnitude lower than the one 
given in Ref. [Q. The discrepancy can be easily understood. We have used a more realistic spherical model, where 
the flux varies as (Eq. (|^), while in Ref. |^ the resonance region was described in terms of a plane Eddington 
atmosphere, where F is constant. In the last case hp^ vanishes and no kick is generated at the lowest order. 
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IV. THE POLYTROPE MODEL 



The conditions at the inner core of the protostar are consistent with a polytrope gas of relativistic nucleons with 
an adiabatic index F = 4/3 In this section we assume that this value of F also holds in the rest of the star. 

This model satisfies the same equations of hydrodynamic equilibrium (|3^) , energy transport (^) and flux conservation 
( ^ ) as the Eddington model, but differs in the equation of state. It leads to a confined atmosphere, where the density 
becomes zero at a radius of the order of a few core radii. 

The equation of state for a polytrope gas of adiabatic index F is given by p^Jl5| , p7[ | 



P ^ Kp^ . 

From Eqs. (^ ) and (^o|), the density profile is determined by the integro-differential equation 

d/'i _ XtRcpI-^ M{r) 



dr 



Mr, 



where Ar — J^^v- 1 ■ extremely good approximation for p is given by the function 



r-i 



Ar ( — - 1 ) m(r) + 1 



(47) 



(48) 



(49) 



with m(r) = /i + (1 — ^) such that m{Rc) = 1 . The radius of the star, defined by the surface where the density 
becomes zero, p{Rs) — 0, determines the value of the p, parameter. The relationship between p and Rg is given by 



R. 



Rr 



Ar (i?s — Rc) Rs 
It is convenient to rewrite Eq. (^ ) as follows 



Rs 



{Rs - Rc) ■ 



r-i 



Rr 



Rr 



a —] +b—+C 



(50) 



(51) 



with a = (1 — /i)Ar, b — {2p — l)Ar, and c = 1 — p\r- 

Once the mass density is known, the temperature profile is determined by the condition of constant luminosity (|34|), 
together with (Eq): 



dT^ 
dr 



9kLc 
2~P 



with p given by Eq. (]49|). The sohition to Eq. (^2|) for T ~ 4/3 can be found exactly. We write it as 



T{T)=Tj\r 



where x(a;) is a polynomial of degree seven in the variable x 



x(-)-x(i) + i 



(52) 



(53) 



3 b 
■)^[x) = c^x + -bc^ x^ + c{ac + b^)x^ + -(6ac + 6^)a;** 

3a, ,9x ba^ r. ^ 
+ — (ac + b^)x^ + —x^ + —x^. 



(54) 



where a, 6, and c are the parameters introduced in Eq. (51). 

The inverse characteristic lengths of the temperature and density at the resonance can now be calculated from Eqs. 
(p3|) and (|4^), respectively 



^A. 



Pr (Tc\ R. 



TrJ Rl ' 

1/3 

2a + & 



Rc\ R 



Rr Rl 



(55) 
(56) 
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To estimate the magnitude of B we use the same values for the core parameters as in the Eddington model. The 
constant K is fixed by the condition K = P/p^ , which gives K = TjrrinPc^^^ = 5.6 x 10"^ MeV^/^. Similarly to the 
previous section, we adopt here Rs = 5.8Rc from a numerical estimation. The radius of the surface of resonance can be 
calculated from Eq. ( ^T]) and the result is Rr = 4.2i?c. Using these values in Eq. ( ^6|) we obtain h^^ = —4 x 10~^R~^ 
and h-^ = -UR^^, which substituted in Eqs. (|l|) and @ yield B ~ 5 x 10^^ in order to have Ak/k ~ 0.01. 
This result is in agreement with the values for the magnitude of the magnetic field calculated in the previous section. 



V. CONCLUSIONS 



In this paper we revisit the resonant neutrino conversion for pulsar kicks. By expressing the momentum asymmetry 
in terms of the logarithmic derivative of the energy flux, we make clear that a kick is produced at order 5/Rr by the 
radial dependence of this flux combined with a deformation of the resonance surface. An important ingredient for this 
result is the existence of a component of the neutrino flux normal to the resonance surface, which acts an effective 
tau neutrinosphere. This is valid even though the neutrino luminosities are controlled by the core emission. In the 
particular case of a plane atmosphere with constant flux, hj,^ vanishes and there is no kick effect to lowest order. 
However, with a more realistic spherical geometry for the atmosphere this is not true anymore. 

To estimate the neutrino flux anisotropy, we consider two simple self-consistent models for the stellar atmosphere: 
the spherical Eddington model and a model for a star composed of a polytropic neutron gas. Both models take into 
account the energy flux conservation and give reasonable profiles for temperature, density, and pressure. The results 
clearly show that the main effect is due to the geometrical deformation of the neutrinosphere surface, and not to the 
alteration of a conserved plane-parallel flux by a response of the global stellar structure |^ . This last contribution 
is much weaker than the geometrical one. It should be noticed that since sterile neutrinos do not interact with 
the medium, the in-going Pg produced at the resonance surface would not be absorved. Thus, in this case, in the 
approximation we use, the resonant conversion mechanism is not able to generate a kick. 

For typical values of the protoneutron star parameters, in both models the magnetic field required to generate an 
appropriate kick is of order of i? ~ 10^® G. A more reliable quantitative evaluation of the effect would require a 
detailed model and much more involved numerical calculation. Nevertheless, our simplified discussion indicates that 
the neutrino oscillation mechanisms can not be discarded as a possible explanation for the pulsar velocities. Here, 
we adopt the usual mechanism of oscillations between massive active neutrinos, which requires a large tau neutrino 
mass of the order of 100 eV . However, the same analysis can be performed in other scenarios in agreement with the 
present boundaries on the neutrino properties, such as oscillations induced by a violation of the equivalence principle. 
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APPENDIX: ENERGY FLUX ON A DISTORTED RESONANCE SURFACE 



Here we analyze the energy flux through the distorted surface of resonance, a sphere with its center shifted in 5 
with respect to the center of the star. At a given point on this surface characterized by a polar angle -d, referred to 
the direction of the magnetic field, the normal forms an angle a with respect to the radial direction, given by 



sin a — — sin d 
R 



(57) 



The expression for the neutrino distribution in the diffusive approximation takes the form 



U{r,k,n) = f:^iv,k) 



1 



Kpk'^ 



eq 



dr 



COS (3 , 



(58) 



where /3 is the angle between k — ktt and the radial direction f . We assume null chemical potential. Introducing 
local coordinates given by an azimuthal angle x on the tangent plane (with % • fi = 0), and a polar angle 4> respect to 
the interior normal (— n), /? can be expressed as 
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COS /? = sin X sin 4> sin a + cos (j) cos a 



(59) 



At a given point of the resonance surface the produced tau neutrinos have the same isotropic distribution as 
the parent electron neutrinos. The in-going tau neutrinos that fall into their own neutrinosphere are absorbed and 
thermalized. The remaining tau neutrinos escape from the star and contribute to the kick. In the following we will 
assume for simplicity that all the in-going tau neutrinos are absorbed. Thus, from Eq. restricting the integration 
to k • n > 0, we obtain that the flux has the normal component {Fn n) and the tangential one {F^ d) given by 



Fa = —<jT^ 
" 64 

1 1 

~ T2~p 



1 1 
144 up 



dr 



cos a , 



dr 



sma , 



(60) 
(61) 



where all quantities are evaluated at the point of the resonance surface defined by the angle Alternatively, we can 
describe this flux in terms of a normal component and a radial one, which are 



7 

F- - —nT'^ 
^" " 64 

1 1 



F, 



144 up 



dT^ 



dr 



(62) 
(63) 



respectively. The normal flux comes from the contribution of /^^^ to the neutrino distribution function, while the 
radial flux comes from the term which depends on its gradient. We observe that the radial flux above the limit surface 
is one half of the flux immediately below that surface. 



dr 



(64) 



which is purely diffusive. Because cosa = 1 + 0{6^ /R^), the equation V • F = 0, integrated in an infinitesimally thin 
volume just above the limit surface, implies that F^ is equal to Ff. Hence, for a given -d we have 



Ff{d) = Ff^id) = ^F,{Rr + 6cosd) ~ ^F,{Rr) (l + hp^S cosd) , 



(65) 



where hj.'^ = -2R-^ . 
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